Draft version February 5, 2008 

Preprint typeset using I^TgX style emulateapj v. 10/09/06 



O 

o 

cn 



THE HARD X-RAY TAILS IN NEUTRON STAR LOW MASS X-RAY BINARIES: Be/?/?o5AX OBSERVATIONS AND 
POSSIBLE THEORETICAL EXPLANATION OF THE GX \l+2 CASE 

Ruben Farinelli', Lev Titarchuk''^-^ & Filippo Frontera' '* 
Draft version February 5, 2008 

ABSTRACT 

We report results of a new spectral analysis of two Be/:i/7o&4X observations of the Z source GX 17-1-2. In one 
of the two observations the source exhibits a powerlaw-like hard (> 30 keV) X-ray tail which was described 
in a previous work by a hybrid Comptonization model. Recent high-energy observations with INTEGRAL of 
a sample of Low Mass X-Ray Binaries including both Z and atoll classes have shown that bulk (dynamical) 
Comptonization of soft photons can be a possible alternative mechanism for producing hard X-ray tails in such 
systems. We start from the INTEGRAL results and we exploit the broad-band capability of BeppoSAX to better 
investigate the physical processes at work. We use GX \l+2 as a representative case. Moreover, we suggest 
that weakening (or disappearance) of the hard X-ray tail can be explained by increasing radiation pressure 
originated at the surface of the neutron star (NS). As a result the high radiation pressure stops the bulk inflow 
and consequently this radiation feedback of the NS surface leads to quenching the bulk Comptonization. 

Subject headings: stars: individual: GX \l+2 — stars: neutron — X-rays: binaries — accretion, accretion 
disks 



> 

\o 

(N 
O 
l> 

o 

Oh 
I 

o 



1. INTRODUCTION 

One of the most interesting features of the high-energy 
X-ray Astronomy since the BeppoSAX and RXTE era is the 
observation of a transient powerlaw (PL)-like emission above 
30 keV in low mass X-ray binaries (LMXBs) hosting a 
weakly magnetized neutron star (NS). This component, which 
appears over the persistent continuum composed by a thermal 
Comptonization (TC) plus a soft blackbody (BB)-like emis- 
sion, up to now has been observed in the six X-ray sources 
belonging to the so-called Z class (GX 17+2, Di Salvo et al. 
2000, hereafter DSOO, Farinelli et al. 2005, hereafter F05; 
GX 349+2, Di Salvo et al. 2001; GX 340+0, Lavagetto et 
al. 2004; ScoX-1, D'Amico et al. 2001; GX 5-1, Asai et 
al. 1994, Paizis et al. 2005; Cyg X-2, Frontera et al. 1998, 
Di Salvo et al. 2002), in the peculiar source Cir X-1 (laria 
et al. 2001, 2002) and, very recently by INTEGRAL (Paizis 
et al. 2006, hereafter P06), in the source GX 13+1. It is 
worth noting, although GX 13+1 is classified as atoll source 
by Hasinger & van der Klis (1989), shows, in fact, spectral 
properties similar to Z sources. 

The hard X-ray emission seems to be correlated with lower 
M and radio-loud states (see, e.g., Penninx et al. 1988 for 
Sco X-1 and GX 17+2), but its physical origin is still debated. 

The thermal Comptonization of soft photons in very high- 
temperature plasma does not seem to apply because it is very 
difficult to explain where and how such a hot plasma could 
be present in a source where a spectral signature of a low- 
temperature (~ 3-5 keV) plasma is definitely observed (see, 
e.g., review by Di Salvo & Stella 2002). 

The first attempt to describe the hard X-ray tail in terms of 
some physical self-consistent model was performed by F05 
who fitted the BeppoSAX broad-band spectrum of GX 17+2 
using a hybrid Comptonization model (Coppi 1999) plus a 
soft ~ 0.6 keV BB-component. The presence or absence of 
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the hard tail was explained by decreasing the power (supplied 
by some unspecified mechanism) injected to accelerate some 
fraction of electrons over their Maxwellian energy distribu- 
tion. 

Titarchuk, Lapidus & MusHmov (1998), hereafter TLM98, 
showed that the Keplerian disk flow should adjust to the sub- 
Keplerian rotation of the central object, NS or Black Hole 
(BH). For a given viscosity in the disk the location of the ad- 
justment radius Rgdj depends on the mass accretion rate only. 
Thus in the transition layer (TL) between Rf^s and Rgdj the 
rotational velocity is sub-Keplerian. Moreover Titarchuk, & 
Farinelli (2007 in preparation), hereafter TF07, calculate the 
radial velocity vr in the TL and they find that vr is very close 
to the sound speed vs in the zone between Rgdj and some ra- 
dius Rff below which the accretion onto NS occurs in almost 
free-fall manner. 

Thus TF07 unambiguously show the bulk inflow must be 
considered a fundamental process in accretion into the com- 
pact objects. Titarchuk, Mastichiadis & Kylafis (1996), here- 
after TMK96, studied the effect of the bulk motion Comp- 
tonization (BMC) of soft photons by energetic electrons in a 
converging flow into NS. They found that the BMC emergent 
spectrum has a specific hard tail which can be a possible ex- 
planation for the hard X-ray emission in NS sources. 

It is worthwhile to emphasize that this BMC process was 
not considered by F05 a possible mechanism to produce the 
observed hard tail of X-ray spectra from NS sources. 

However P06 analyzed high-energy (> 20 keV) observa- 
tions of NS systems performed by INTEGRAL and they con- 
firmed that BMC indeed can explain this hard X-ray emission. 

P06 supported this spectral evolution picture of NS sources 
using a sample of twelve bright low-mass X-ray binaries host- 
ing a NS. Their sample comprises the six Galactic Z sources 
and six atoll sources, four of which are bright GX bulge 
sources while two (H 1705^40 and H 1608-552) are weaker 
in the 2-10 keV range. 

Comparing their results with those obtained by Falanga et 
al. (2006) from GX 354-0 (4U 1728-34), P06 identified four 
main spectral states for NS LMXBs (see Fig. 4 in P06): 
low/hard state (GX 354-0), hard/PL state (H 1705-440 and 
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H 1608-552), intermediate state (e.g., GX 5-1) and very soft 
state (e.g., GX 3+1). Above 20 keV, the low/hard spectra 
are described by TC with plasma temperature kT^ ~ 30 keV 
and optical depth t ~ 1; the hard/PL state spectra, as sug- 
gested by the name, show simple PL-like emission; the in- 
termediate state spectra are characterized by the presence of 
a TC-component with spectral parameters similar to those of 
the soft state plus additional PL-like emission (above ~ 30 
keV); and the soft state spectra are described by a single TC- 
component with low-^T^e and high-r. 

We remind the reader that P06 used the high-energy (> 20 
keV) INTEGRAL observations of these NS sources, so they 
could not reveal the soft BB-component observed in NS 
LMXBs (see review by Barret 2001). Using broad-band data, 
the intermediate states of NS (see, e.g., GX 17+2 and Cyg 
X-2), being similar to high-soft states of BHCs, clearly show 
the presence of a soft BB-like component plus a PL emis- 
sion at high energies, but with an additional TC feature that 
is absent in BHC spectra. This is because the presence of a 
solid surface leads to a further channel of emission. Thus, this 
similarity of spectral appearances in NSs and BHCs suggests 
a similarity in the physical processes occurring in these two 
classes of systems. 

In this work, using the excellent broad-band observing ca- 
pability of BeppoSAX (Boella et al. 1997a), we present the 
0.4-120 keV energy spectrum of the Z source GX 17+2, in 
order to shed light on the source geometry. We find that the 
BMC model can really explain the origin of the high-energy 
emission in NS LMXBs. 

2. THE BULK COMPTONIZATION MODEL 

The TC theory along with the dynamical (bulk motion) 
Comptonization has been first suggested by Blandford & 
Payne (1981, hereafter BP81), and further developed in de- 
tail by Titarchuk, Mastichiadis & Kylafis (1997, hereafter 
TMK97). BP81 derived the diffusion Fokker-Planck equa- 
tion using the sub-relativistic kinetic equation [see Eq. [15] in 
BP81 and Eq. [14] in TMK97]. The relative efficiency of TC 
effect with respect to BMC is determined as (see TMK97) 



< hE,h > 1 



(1) 



where < Ais,/, > and < AEi, > represent the average energy 
photon exchange due to TC and BMC respectively. The 6 
parameter is defined as (see TMK96 and TMK97) 
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where = kTs/nieC^ is a dimensionless plasma temperature 
and TiQ = kTJ(lO keV), while i = L/LEdd and rit = M/MEdd 
are a dimensionless luminosity emitted by the NS and mass 
accretion rate calculated in Eddington units (Mndd = ^Edd/c^) 
respectively. When 6 = 0, diffusion equation (14) of TMK97 
is reduced to the TC difiiision equation (see, e.g., Titarchuk 
1994). For our purposes it is just sufficient to point out that 
the photon occupation number of the emergent spectrum can 
be expressed as 

F(x) = ^ I Ik{x,xo)sixo)dxo = ^ Ik{x,xo) * sixo), (3) 



<:=1 



k=l 



where x = E/kT^, xq = E/kT^b, kT^, kT^b are temperatures of 
plasma and of seed photons, respectively and * is the convo- 
lution operator. 



TMK97 showed that /i (x, xq) oc x " for A; = 1 and x > xoup 
to a rollover energy of order of Ec ~ nieC^/m, while 4(x, xq) oc 

/>oo 

5(x-;ico) for ^ > 2, so that /(t*.s oc 6(x—X(y)s(xo)dxo ^ s(x). 
In other words, only the first term (k = 1) in the series of 
Eq. ||3| significantly contributes to the upscattering (Comp- 
tonization) part of the emergent spectrum. The terms related 
to ^ > 2 represent the photons which escape the plasma cloud 
with insignificant energy-change (see also Fig. 3 in TMK97) 
and thus they retain the shape the original seed photon spec- 
trum. 

The relative importance of the Comptonization term /i * 
s with respect to that of the seed photons depends on their 
spatial distribution within the cloud (see Fig. 4 in TMK97). 
In the BMC model for xspec developed by TMK97, the seed 
photon spectrum s{x()) of Eq. ||3] is assumed to be a BB-like 
spectrum, while the relative importance of the efficiently up- 
scattered photons, the term /] * BB of series (O with respect 
to the direct BB component is expressed using a weighting 
factor A/(A + 1). In fact, the xspecBMC model reports log(A). 
Thus the analytical formula (back-reverting from x to E) for 
the total spectrum is 



F(E) 



N 



A+ 1 



(BB+A Ii* BB). 



(4) 



The free parameters for xspec are temperature of the seed pho- 
tons kTbb, spectral index a of the term /i, a logarithm of the 
weighting factor log(A) and the normalization Cn- The spec- 
tral index a is related to the efficiency of the Comptoniza- 
tion processed (both thermal and bulk). As higher cc-value, 
the Comptonization efficiency is less and thus the system is 
closer to thermal equilibrium, so that the emergent spectrum 
has a BB-like shape. The weighting factor A is chosen in such 
a way that for A — > the model is reduced to the standard BB 
model of xspec 

3. RESULTS 

Using XSPEC vll.3 we made a spectral analysis of a set of 
two BeppoSAX observations of the Z source GX 17+2 per- 
formed in 1997 (on April 3 and 21, respectively). A detailed 
description of the data processing and analysis can be found in 
F05. We just point out that in both cases the source was in the 
horizontal branch (HB) of its colour-intensity diagram (HID). 
As during the first observation the source stayed in a clus- 
tered region of the HB, whereas in ffie second observations 
the source evolved along the branch, from right to left (see 
Fig. 2 in F05). The HB region of the second observation was 
divided into three parts, for which separated spectral analysis 
was carried out. Given that there was only a very slight varia- 
tion of the spectral parameters in these regions, F05 reported 
their results in terms of the time-averaged spectrum. 

In the first observation F05 found a systematic excess in 
the residuals above 30 keV by fitting the persistent continuum 
with a photoelectrically-absorbed (wabs in xspec) BB plus the 
TC model of Titarchuk (1994, comptt in xspec). This excess 
was phenomenologically described by a PL with photon index 
r ~ 2.8. DSOO obtained the same F value for a different set of 
BeppoSAX observations, when the source was in the left part 
of the HB. 

Given that the output shape of the physical BMC spectral 
model consists of a direct (non up-scattered) BB-like com- 
ponent along with a PL shape at high energies, our new fit- 
ting model was comptt plus bmc plus a gaussian emission line 
around 6.7 keV, which was previously observed both by DSOO 
andF05. 
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To vary all of the model parameters, we find that the data 
can be presented by two solutions in the framework of the 
same model. One solution is with kT\,\, < kT^ (Case A) and 
another one is with kT^h > kT^ (Case B). We remind a reader 
that kTbi, and kT^ are seed photon temperatures for bmc and 
coMPTT models respectively. 

We note that this dichotomy was already reported by F05 
for the BB plus comptt model. In both cases we obtain 
;^'^/dof= 151/144. It implies that we cannot statistically dis- 
tinguish between these two solutions. In Case A we find 
/trbb = 0.62 + 0.02 keV, kT^ = 1.14 ± 0.05 keV, log(A) = 
0.05 (> -0.66), and a = 2.1+°^^ while for Case B the best-fit 
parameter values are reported in Table[T](see the first column). 
It is worth noting that both the comptt and Gaussian best-fit pa- 
rameters are almost the same within error bars for these two 
solutions. 

In the second observation, the source was detected {S IN > 
3) up to 50 keV. A simple BB plus comptt plus Gaussian model 
provides ;(^2/dof=192/142 and ;if2/dof=197/142 for Case A 
{kTbh < kTv/) and B {kT\,b > kT„) respectively. As already 
discussed in F05, the reason for these non excellent val- 
ues can be attributed to averaging different (even though very 
close each other) spectral states. Indeed, no systematic devi- 
ations are observed in the residuals, and the last PDS point 
(40+ 10 keV) is at ~2.5 cr level above the model. As a con- 
sequence of this, the bmc parameters a and log(A) remains 
largely undetermined if they are left free in the fit. We thus 
fixed a to 2. 1 and 2.6 for Cases A and B respectively. In both 
Cases the fit quality only slightly improves, the corresponding 
X^/dof decreases to value of 187/140. 

A comparison between the two observations, shows that the 
major change is related to log(A), whose value decreases from 
~ 0.73 in the first observation to 0.33 in the second obser- 
vation. This is consistent with the fact that the PL-like emis- 
sion is significantly detected only in the first observation (see 
Table 1). 

Below we present a physical reasoning why the best-fit pa- 
rameters of Case B are more preferable than that for Case A. 
In Fig. [r|(top panel) we present the deconvolved EF(E) spec- 
tra of the two observations for Case B to illustrate this eff'ect 
of the transient hard X-ray emission. 



TABLE 1 

Best-fit parameters of the multi-component model wabs(bmc + comptt + 
GAUSSIAN). Errors are computed at 90% confidence level for a single 

PARAMETER. 



Parameter 


Obs. 1 


Obs. 2 




1 qi+0.07 


1 89+"-"'' 

'■°'-0.08 


BMC 

kTu (keV) 
log(A) 
a 


1 33+0.14 

-0.10 

0.73 (> -0.24) 

2 64+" ''' 


1 31+0.07 
— ^ -0.08 

-0 33+°-29 

■ -0.39 

[2.64] 


COMPTT 

kiw (KeV) 

Kl ^ \I^C V f 
T 


0.54 

3 39+0.07 
■^■■^ -0.06 

12 5+" '' 
''^■-^-04 


fi c^+0.03 

■ -0.03 
3 1 1 +0.09 

■ -0.08 

13 1+"* 

10.1 Q J 


GAUSSIAN 

El (keV) 
o-i (keV) 

A* 

EW, (eV) 


6 69+" °5 

"■^'-0.09 

38+« 


"■' -0.09 

24+"-" 

"■^^-0.12 

^■i^li 

35!" 




0.16 


0.21 



L"/ 1.46 1.47 

X^/doi 154/145 185/141 

" In units of 10^^ cm"^. 

Total photons in the line in units of 10"^^ cm"^ s"' . 

Extrapolated in the energy range 0.1-200 keV. 

In units of 10^* erg s"' assuming a distance of 7.5 kpc 
(Penninx et al. 1988). 



4. DISCUSSION AND CONCLUSIONS 

In this Paper we show the results of a new X-ray spectral 
analysis, using BeppoSAX broad-band data, of the Z source 
GX 17+2 during two states in which a PL-like hard X-ray tail 
was present and absent, respectively. We describe this hard 
tail in the framework of Comptonization of soft BB-like seed 
photons in the converging flow into NS. 

Comparison of our best-fit parameters with those reported 
in P06 reveals some significant differences: Restricting 
energy-band to above 20 keV allows one to see only the ther- 
mal (BB-like or Wien-like) bump plus the additional PL-like 
X-ray tail in the spectrum (see, e.g.. Fig. 5 in P06). This ther- 
mal bump, lacking low-energy information, is interpreted by 
BMC as the direct seed photon component, namely the com- 
ponent which consists of the photons subjected insignificant 
energy change when they scatter in the ambient medium. This 
is confirmed by the fact that the color temperature ^T^bb of this 
component about ~ 2.7 keV (see Table 2 of P06), is not far 
away from the electron temperature of comptt found by us 
(~ 3.3 keV, see Table[B. 

However using broad-band data it is possible to see an addi- 
tional direct seed thermal emission located at lower energies 
whose thermal Comptonized part dominates in the 20-40 keV 
region (see Fig. [T] top panel). One unavoidable consequence 
of these considerations is also the different values of the illu- 
minating factor A [or log(A)]. 

The broad-band analysis allows us to elaborate a more de- 
tailed picture of the source properties. As shown in Sect. [3] 
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we can identify two regions of seed photon temperatures, for 
which the temperature values differ from each other by a fac- 
tor about two. In Case A the hotter (~ 1.2 keV) photons are a 
source for thermal Comptonization, while the cooler ones (~ 
0.6 keV) are up-scattered to form the PL emission at high en- 
ergies. In Case B the opposite situation occurs. Given that we 
cannot statistically distinguish one model from another, the 
models give the same x'- value, only considerations based on 
the physical characteristics of the emitting regions can help us 
to discriminate among them. 

Several theoretical studies (e.g., TLM98; Popham & Sun- 
yaev 2001, hereafter PSOl) have shown that when the geo- 
metrically thin accretion disk approaches the NS, the differ- 
ence in the angular velocity of the disk and the more slowly 
spinning accreting star gives rise to a transition region, usu- 
ally called the boundary layer (BL) in PSOl or transition layer 
(TL) in TLM98. PSOl argue that when M increases, the BL 
extends radially (reaching 2-3 /?ivs), significantly increasing 
its contribution to the total emitted luminosity (even to 70%). 
Moreover PSOl show that inside the BL, the plasma proper- 
ties (e.g., temperature, density, pressure, radial velocity) are 
not constant but vary as a function of the distance from the 
NS surface. 

TF07 demonstrate that in the TL, the radial velocity profile 
as a function of radius is almost constant from the adjustment 
radius Radj (defined as the outer TL radius) to some lower 
radius Rff where the constant velocity flow is followed by 
the free-fall accretion flow, affected by radiation pressure (see 
also Introduction section). TF07 obtain this radial velocity 
profile as a solution of the radial momentum equation [ see 
Eq. [5] in PSOl] for which they use a solution of the angu- 
lar momentum equation found by TLM98 (see also Titarchuk 
& Osherovich 1999). The free-fall region of the TL is likely 
a place where the bulk motion Comptonizarion (BMC) spec- 
trum is formed. 

The free-fall region is the innermost part of the TL, and con- 
sequently the seed photons subjected to BMC mainly come 
from the NS surface. The temperature of the NS photons 
kT\jQ is higher than that of the disk soft photons kT^f, be- 
cause the effective area of NS surface ~ 3/?s (^s is the 
Schwarzchild radius) is more compact than that of the Keple- 
rian disk ~ (10 - \5)R^ (see Shrader & Titarchuk 1999). This 
is a physical (not statistical) reason why we choose, among 
the two solutions, the one for which kT\jo > kT^. In Fig. [T] 
(bottom panel) we provide a schematic view of the possible 
source accretion geometry. 

The conservation of energy provides a free-fall radial ve- 
locity given by 



vg^c[(l-F/FEiA)(Rs/R)] 



1/2 



where F is the local flux and 



GMcmnnijii 
R-criric 



(5) 



(6) 



is the local Eddington flux. In Eq. IS), cr-j- is the Thomson 
cross-section, M is the NS mass, c is the speed of light, n,- 
and M, are the ion density number and molecular weight, re- 
spectively, while «e is the electron density. As it follows from 
Eqs. ||2l and ||5| Vff oc ^ . When the local flux approaches 
the local Eddington flux, the factor (1 - f/^Edd)'^^ goes to 
zero and consequently Vff goes to zero too. Thus the eff'ect of 
the disappearance of the BMC hard tail {5 — > 0) can be natu- 
rally explained as a result of the feedback from the NS surface 



when the NS emergent (local) radiation flux is very close to 
Eddington. 

In Table [1] one can see that, while the total flux remains al- 
most the same, the contribution of the seed BB-component 
(suggested to come from the NS) increases from the first ob- 
servation (Fbb/-Ftot = 0.16) to the second one (Fbb/-Ftot - 
0.21). Given that the BB-flux flux is proportional to the to- 
tal M, we can conclude that M, is a quantity which regulates 
appearance and disappearance of the BMC hard tail in X-ray 
spectra from NS LMXBs. 

If the firm surface is absent, like in BHs, then for a given 
electron temperature of the converging bulk flow the index F 
decreases with M to some value and then it saturates because 
of the photon trapping effect in the converging flow (see more 
details in Titarchuk & Zannias 1998). But in the NS case, as 
we explain above, there is a feedback from the NS surface 
when M is very high. The high radiation flux emerged from 
the NS surface quenches the bulk free-fall flow. However, we 
may not exclude that the hard tail, instead of disappearing, 
just decreases its intensity, going below the instrument detec- 
tion threshold, because of a lower fraction of seed (NS) pho- 
tons intercepted by the converging flow {geometrical eff'ect). 
The amount of weakly upscattered seed photons leaving the 
flow depends on their spatial distribution within the flow it- 
self (TMK97). The BB-like direct component is stronger for 
more uniform spatial distribution of seed photons than that 
mainly concentrated at the bottom. 

The changes of the observed contributions coming from the 
BB-like and the PL component may be equally provided ei- 
ther by a change of the seed photon spatial distribution (a 
more uniform distribution would decrease the amount of effi- 
ciently up-scattered photons, spatial effect) or by a lower sub- 
stended angle of the flow as seen from the NS (geometrical 
effect). Another important point of our discussion concerns 
the apparent mass accretion rate M. Very often, changes in 
M in a given sources are inferred by measuring changing on 
its apparent luminosity, using the relation M= F^4nD^/(ric^). 
However, this approach may lead to controversy. In fact, Di 
Salvo et al. (2002) and DSOO observed the hard X-ray tail 
in Cyg X-2 and GX 17+2 when the sources were in the HB 
of their HIDs. The apparent 0.1-100 keV luminosity L„pp 
evolves during the source motion along Z-track from the HB 
to the normal branch (NB). Lapp increases in Cyg X-2 and it 
decreases in GX 17+2 with S- (where the length parameter 
5 J along the Z-track increases from HB to NB). 

Moreover the NS sources belonging to the high-luminosity 
atofl class GX 3+1, GX 9+1 and GX 9+9 have stable X- 
ray spectra very close to Z sources but, up to now, no evi- 
dence of PL-like hard X-ray emission has been provided for 
them. This is in constrast with what expected, given that, de- 
spite their spectral similarity, these sources are less luminous 
than Z sources. Additionally, our data analysis shows that the 
hard X-ray behaviour changes at the same apparent luminos- 
ity level (see Table [U. All these observable quantities, show 
that mapping of M changes either within the same source, or 
among different sources must be treated very carefully. 

Thus we can suggest that the NS surface feedback on the 
BMC tail is not sensitive to the absolute (disk plus TL) M 
value, but it mostly depends on the local radiation flux related 
to the Mtl impinging on the converging flow. This value may 
significantly differ from the simple L/(4nR~) estimate (see 
also PSOl). 

From the results of the multi-component fit, we are also able 
to map the energy budget of the system. We reveal that a frac- 
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tion of the low-temperature (0.55 keV) BB radiation in the to- 
tal flux, which is presumably the disk contribution to the total 
flux FdisklF,ot, is only ~ 34%. The larger fraction of the grav- 
itational energy of the accreting material ~ 66% is released in 
the transition layer and at the NS surface. We infer FaisklFtot 
using a fraction of the TC-component (comptt) in the total 
flux Fcomp/F,oi ~ 78% and a value of the Compton enhance- 
ment factor £comp ~ 2.3. Note the values of Fcomp/^wr and 
^comp are almost the same for the two observations. We re- 
mind a reader that Scomp is defined as a ratio of the flux of 
the Comptonized emerging spectrum Fcomp with respect to 
the soft photon injected flux Fsofi (see Sunyaev & Titarchuk 
1980). We calculate Ecomp using coMPTT-model. Namely, the 
model flux with the best-fit coMPTT-parameters gives us Fcomp 
whereas the model flux with the best-fit kT^, kT^, normaliza- 
tion and very small optical depth t provides us Fsoft^- 

We want also to emphasize that our preference to explain 
the transient hard X-ray tails in NS LMXBs in the frame- 
work of the BMC phenomenon is not related to the fact that 
this model better describes the data statistically. The hybrid 

^ Given that the average number of scatterings for photons is oc t^, it is 
just sufficient to set t <*: 1 in comftt and obtain Fs„f,. 



Comptonization model by Coppi (1999) first used by F05 to 
explain the GX \l+2 BeppoSAX broad-band spectrum was 
statistically acceptable too. 

The main reason for our choice is that the BMC process is 
just related to the dynamics of the accretion flow, based on the 
exact solution of dynamical equations, and does not require 
any particular tuning and any particular energetic configu- 
ration of the Comptonizing plasma (e.g. thermal plus non- 
thermal distribution). We thus cannot a-priori exclude that 
non-thermal electron distributions are present in such a sys- 
tems (for instance in possible emitting jets), but the key point 
is that they are not required to explain our data. 

In BH sources, the BMC phenomenon are likely even more 
prominent than in NSs (because of the absence of a solid sur- 
face and thus the absence of a strong radiation feedback effect 
on the accreting flow) and their effect can well explain the 
transition from hard/low to high/soft states in these systems. 

We acknowledge the referee for his/her important questions 
and for the constructive suggestions which improved the pa- 
per quality. 
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Fig. 1. — Top: Unabsorbed EF(E) spectra and best-fit model wabs(bmc + 
COMPTT + GAUSSIAN) for the first {upper panel) and the second {lower panel) 
observation of GX 17+2 performed by BeppoSAX. The residuals (in units 
of lo") between the data and the model are also shown below each spec- 
trum. Dift'erent line styles represent the single components. Dotted-dashed: 
COMPTT. Dotted: bmc. Long-dashed: Gaussian. Bottom: Schematic view 
of the accretion geometry of GX 17+2 as inferred from the single spectral 
components observed in its X-ray spectrum (see Top). The decreasing (or 
disappearing) of the PL-like emission at high energies may be due to a de- 
crease of the substended solid angle of the radiation-corrected free-fall bulk 
region as seen from the NS, or to higher radiation pressure. It is also possible 
that both effects lead to the hard tail auenchine. 



